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The dynamics of the magnetic order parameters of ferromagnet/paramagnet/ferromagnet spin 
valves and isolated ferromagnets may be very different. We investigate the role of the nonequilibrium 
spin-current exchange between the ferromagnets in the magnetization precession and switching. We 
find a (low-temperature) critical current bias for a uniform current-induced magnetization excitation 
in spin valves, which unifies and generalizes previous ideas of Slonczewski and Berger. In the 
absence of an applied bias, the effect of the spin transfer can be expressed as magnetic-configuration- 
dependent Gilbert damping. 



PACS numbers: 75.40.Gb,75.70.Cn,76.60.Es,76.50.-|-g 

The giant-magnetoresistancei (GMR) in multilayers of 
metallic ferromagnet {F) and normal-metal {N) films 
controls the current flow by a magnetic-field-induced re- 
orientation of the magnetizations. This effect has found 
applications in, e.g., magnetic-field detectors as read 
heads of mass storage devices. Modern magnetic stor- 
age media and magnetic random-access memories also 
consist of F/ N composites with information stored by 
switchable magnetic configurations. Device performance 
is measured in terms of bit density as well as speed of 
reading and writing information. A thorough under- 
standing of the magnetization dynamics in F/ N hybrid 
structures is therefore desirable. 

Magnetization reversal is usually achieved by magnetic 
fields generated by external electric currents. In small 
structures, however, much energy is wasted in the form 
of stray magnetic fields, which motivates consideration 
of other switching mechanisms. An effect inverse to the 
GMR (i.e., exerting spin torques on the magnetizations 
by an electric current passed through the multilayer) ap- 
pears rather promising for this purpose. The spin torque 
can lead to magnetization switching at a critical electric 
current as has been experimentally confirmed^i^ 

Perpendicular spin valves, i.e., Fg/ N/ Fh trilayer pillar 
structures with layer thicknesses down to a few mono- 
layers and lateral dimensions in the (sub)micron region, 
are ideal to study precession and switching phenomena 
in hybrid systems. When reservoirs are attached on the 
outer sides, these spin valves can be biased by an elec- 
tric current perpendicular to the interface planes. F^ is 
a "soft" ferromagnetic film with a magnetization that 
can change easily, whereas Fh is a "hard" magnetic layer 
whose magnetization is assumed to be stationary. The 
relevant variable is then the time-dependent magnetiza- 
tion of the soft layer. In the following, we will consider 
two scenarios: that of exciting the soft layer by a current 
bias or driving it by an applied rf magnetic field. In the 
former case, the layer F^ can be made stationary by mak- 
ing it much thicker (and therefore more inert for a given 
spin torque) than Fs or by using a resistance anisotropy^ 
while in the latter case (realized as, e.g., an isolated mag- 



netic bilayer), Fh can be pinned by an exchange bias 
or surface magnetic anisotropyr^ For small enough sys- 
tems, the magnetic layers are monodomain ferromagnets 
characterized by two magnetization vectors. Assuming a 
thick enough spacer N, we disregard the static interaction 
between the ferromagnetic layers, but not the dynamic 
coupling^ induced by the spin pumping. 

Slonczewski'^- and Bergei"^ were the first to predict new 
time-dependent effects in spin valves. Both authors have 
realized that a current flowing through a spin valve causes 
a spin transfer through the nonmagnetic spacer, inducing 
spin torques on the ferromagnets. In addition, Berger 
predicted that the two ferromagnets interact via spin 
transfer even in the absence of an applied electric cur- 
rent, resulting in a significant contribution to the Gilbert 
damping of the magnetization dynamics. He further 
demonstrated that a sufficiently large electric current can 
induce coherent spin- wave emission in the ferromagnet, 
an idea which was later supported experimentally^ The 
condition for spin-wave emission^ is similar to the crite- 
rion for the magnetization switching due to Slonczewski^^ 
who treated the Gilbert damping parameter as a phe- 
nomcnological constant. In Ref. a dependence of the 
damping parameter in spin valves on the relative magne- 
tization angle was found. Some of Berger's and Slon- 
czewski's results as well as the underlying theoretical 
models were thus not consistent with each other. In this 
report we offer an alternative theory, which both unifies 
and extends the seminal work of these pioneers. 

Based on the concegt of parametric spin pumping, we 
demonstrated in Ref. llCl that the magnetization motion of 
a ferromagnetic layer is damped by emitting (pumping) 
spins into adjacent conductors. The presence of a sec- 
ond ferromagnet can considerably affect the relaxation 
of the pumped spins, and therefore the magnetization 
dynamics, as discussed below. We combine adiabatic 
spin pumping with magnetoelectronic circuit theorjiiSii^ 
to provide a self-contained framework for spin transfer 
in spin valves. The main results presented here are the 
critical current bias for a low-temperature magnetiza- 
tion instability and the configuration-dependent Gilbert 
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FIG. 1: Schematics of a current-biased spin valve. The sym- 
bols are explained in the text 



damping parameter. In terms of conductance parameters 
accessible to first-principles calculationai^ and combined 
with micromagnetic simulations, the full range of the pre- 
cession and switching dynamics can then be studied in 
principle. 

We consider the system sketched in Fig. ^ The 
Fg/N/Fh trilayer is sandwiched between two normal- 
metal contacts sustaining a charge current bias J. The 
soft layer Fs magnetization m will start moving from its 
equilibrium direction at a critical value Jc (depending on 
the applied magnetic field). Thermal activation facili- 
tates current-induced magnetization switching)& but we 
focus here on the low-temperature regime. The gener- 
alized Landau-Lifshitz-Gilbert equation^- for the mag- 
netization direction m(t) of F^ in the presence of a spin 
current flowing out of F^ reads 
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where 7 is the absolute value of the gyromagnetic ratio 
and Hoft an effective magnetic field which is determined 
by the applied dc magnetic field and the anisotropy en- 
ergy (i.e., shape, crystal, and surface anisotropy). Fg is 
characterized by ao, its intrinsic (dimensionless) Gilbert 
damping constant, M^, its saturation magnetization, d, 
its thickness, and S, its cross section. The spin current 
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consists of the dynamic-exchange current 1°'^'^'' induced 
by the spin pumping and of the current driven by an 
applied current bias. The former has recently been shown 
to be responsible for a dynamic coupling between the 
ferromagnets^ and, according to Ref. |3|, also determines 
the threshold for spin-wave emission. Alternatively, one 
can interpret it as a "dynamic stiffness," which stabi- 
lizes the relative magnetization configuration of the spin 
valve against the torques exerted by I^"''' or an applied 
magnetic field. In high-density metallic systems, the ap- 
plied voltages and spin accumulations are safely smaller 
than the Fermi energies, which means that we are in the 
linear-response regime and both spin currents may be 



calculated independently of each other. Spin pumping in 
the outward direction, i.e., into the external connectors, 
would only increase the intrinsic damping coefficient ao 
by a constant valucjiiS so we disregard it here for simplic- 
ity. 

Let us first consider the spin current IP"'"p pumped 
into the spacer by a time-dependent m(t) in the absence 
of an applied bias, I^"**^ = 0. The transverse component 
of the spin accumulation is absorbed in the ferromagnetic 
layer on the scale of a (transverse) spin-coherence length 



-K/\kp — kp\ in terms of the spin-dependent Fermi wave 

vectors k}p^ (amounting to only a few monolayers for 
transition- metal ferromagnetsiSii^) . When F, is thicker 
than that,^" 

jpuinp 



Here g"^^ — g}^ 
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igl " is the (dimensionless) mixing con- 
ductance of the Fs/ N interface. For simplicity, the con- 
ductance parameters of the two F/ N interfaces are taken 
to be identical in the following. The mixing conductance 
is to a good approximation real- valued, i.e., g]^ <C gj^, at 
least for transition- metal fcrromagnets.^^ jpump ^j-ga^es 
a spin accumulation /Xs^r in N, which induces a back- 
flow spin current Isf^*^ into both ferromagnets i = 1,2. 
According to the circuit theoryjiSiH making use of the 
zero-electric-current condition through the interfaces. 
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Here g'^^ is the (dimensionless) spin-s interface conduc- 
tance, mi — m, ni2 — M, and Afigi = fJ-sN — i^sFi'^i is 
the spin-accumulation difference across the Fi/N inter- 
face. Note that for intermetallic interfaces Sharvin con- 
tributions must be subtracted if conductances are com- 
puted microscopically by scattering theory.'''- The time 
scale of the magnetization dynamics, ~ 10^^^ s, is much 
larger than typical electron dwell times in the metallic 
spacer. Assuming weak spin-flip scattering in A'^, the con- 
servation of spin then implies that 
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The exchange current is given by the difference between 
the pumped current and the backflow: I™"^*^ = IP"™p — 

Tback Tback 

^sl — ^s2 ■ 

The longitudinal component of the spin accumulation 
can penetrate ferromagnets on the scale of the spin- 
diffusion length Asd- In order to find l^i'^^, we solve the 
diffusion equation for the (longitudinal) spin transport 
in the ferromagnets, assuming that the spin current van- 
ishes on the outer boundaries of Fg and Fh ■ It is shown 
that the longitudinal spin-accumulation flow into a fer- 
romagnetic slab of thickness d is then governed by an 
effective conductance g* defined by 
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in terms of.gsd = (^/e^)('S'/Asd)(2(TV''')/(cr^ +(7^'), where 
(7* is the spin-s conductivity of the ferromagnetic bulk, 
so that the backflow current, Eq. I^J, can be written as 
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(7) 

f/* — > when d <^ Asd, i-e., when the spin-flip relax- 
ation vanishes, or when the ferromagnet is halfmetallic, 
so that it completely blocks the longitudinal spin flow 
due to charge conservation. The parameter = {g^^ — 
g*)/{g^^ + g*) characterizes the asymmetry of the ab- 
sorption of transverse vs longitudinal spin currents. Let 
us estimate typical values of v for sputtered Co/Cu and 
Py/Cu hybrids at low temperatures, taking d = 5 nm. 
The main difference between the two combinations is the 
spin-diffusion length in the ferromagnets: Co has a rel- 
atively long Asd ~ 60 nm, while Asd ~ 5 nm is very 
short in Py.^i Using known values for spin-dependent 
conductivities we thus find gsdS~^ « 2.7 nm~^ for Co 



2c,TT£,ii/(c,TT + £,ii)5'- 



20 nm" 



and 16 nm^ for Py. 

for Co/Cu interfacesi^ and we expect the value for Py/Cu 
to be similar. Finally, taking g^^S~^ « 28 nm~^ for the 
Co/Cu interfacei^ and 15 nm~^ for Py/Cuji^ we find 
1/ PS 0.98 for Co/Cu and ly « 0.33 for Py/Cu. 

With I' the same for both layers, we solve for the spin 
accumulation in the normal metal, fisN, in the absence 
of applied current, J = 0, by using the spin conservation, 
Eq. 10, and Eq. O for the backflow in terms of HsN- We 
then finally arrive at 
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M — i^m cos 6 
1 — v"^ cos^ 9 



(8) 



Semiclassically, this equation can be understood as a mul- 
tiple scattering of spin current between the interfaces at 
which the longitudinal part is reflected with probability 
P (X and the transverse component with P oc \ — v. 
We have taken the spacer to be ballistic, so that is 
uniform. Otherwise, the exchange current will be sup- 
pressed by diffuse scattering in the interlayer A^. It is 
straightforward to extend our theory to take this into ac- 
count by, e.g., solving the spin-diffusion equation in the 
spacer and using the same boundary conditions, Eqs. (O 
and Q, as above. 

The magnetization dynamics (in the absence of an ap- 
plied bias) is determined by substituting I™'^*^ into the 
LLC equation, which thus has a damping term that 
cannot be modeled by a constant effective Gilbert pa- 
rameter. We now analyze the configuration dependence 
of the damping in more detail, which is experimentally 
accessible by the FMR line-width broadening at high 
rf intensities^^ (and therefore finite "precession cones" ) . 
For m precessing around M, 
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The angular dependence of the additional Gilbert damp- 
ing parameter due to the exchange spin current then 




FIG. 2: Solid lines are our prediction for the precession-cone 
angle dependence of the Gilbert damping parameter [Eq. dlj] 
and the dotted hnes are Berger's [Eq. llU l. The lower (solid) 
line is representative for Co, while the upper for Py, assuming 
thickness of 5 nm. We expect Fe and Ni to be characterized 
by the two lower (solid) lines. 



reads 



a\e)/a'{Q) = 1 - i/sin^ 61/(1 
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where a'(0) = 'jfigl^ / (SirMsSd) is the damping enhance- 
ment in a coUinear configuration. Interestingly, this re- 
sult bares similarity with Berger's^ 



a'{e)/a'{0) = l/[l + s(l -cosf 



(11) 



where s oc Tsf , a characteristic spin-fiip time, and his form 
of a'{0) is similar as welli^ Expressions ifTUIl and l(TT|) 
are compared in Fig. |21 they have the same small-angle 
asymptotics, but are qualitatively different at large an- 
gles (in particular, the latter is symmetric with respect 
to cos6' while the former is not). At small 0, we can 
rewrite Eq. H10|l to exactly reproduce Eq. (|ll|l after iden- 
tifying s = 1^/(1 — I/). As mentioned above, v is close 
to 0.98 for cobalt and s should be of the order of 100 
(s — 333 is found in Ref. 3 for Co/Cu with Co 1.5 nm 
thick, which remarkably would be quite similar to our 
estimate for this thickness), so that the lower solid line 
in Fig. 121 represents the damping for Co. The "preces- 
sional stiffness" is thus significantly reduced for angles 
which only slightly deviate from the coUinear configura- 
tions (we expect this conclusion to be also true for Fe 
and Ni). Modeling of the magnetization dynamics with 
a constant damping parameter is thus not allowed for 
sufficiently thin magnetic layers. For permalloy, on the 
other hand, the precessional stiffness is expected to re- 
main significant for all angles, see the upper solid line 
in Fig. 121 This implies that the magnetization reversal 
has higher energy-dissipation power, but can occur faster 
than in cobalt in field-induced switching. If m moves 
away from M, i.e., only the relative angle 9 changes, 
then IP"'"P ± M and Eq. © reduces to If^ = IP"'"p/2. 
The "tilting stiffness" has thus an angle-independent en- 
hancement with respect to the intrinsic Gilbert damping. 
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which is given by the same expression as a'(0), i.e., the 
damping in a coUinear configuration. 

Introducing an apphed current bias gives an additional 
control over the magnetization dynamics. When the con- 
ductance parameters of the spin valve are symmetric, the 
bias-induced spin transfer I^'^'' is coplanar with the mag- 
netization directions, I^''"'' = 4^'''''(nH- M)/(2cos6'/2), 
9 being the angle between m and M. This is clear 
after expanding the spin current as I^"*''(m, M) = 
/ii(cos0)m-|-/22(cos6')M + /i2(cos6')m x M and noting 
that I^'^''(m, M) = I^'^"(M, m), which imphes /n = /aa 
and /i2 = 0. The electric current corresponding to a 
given spin-current bias depends on 9 and can be calcu- 
lated readily by circuit theory»iSiii 

Eqs. QJ, (O, 1® and the form of the bias current com- 
pletely determine the dynamics of m(t). The exchange 
induced by the spin pumping causes relaxation toward 
an equilibrium configuration, while the bias current can 
either relax or excite a perturbation from an equilibrium, 
depending on the sign of J. In the process of, e.g., switch- 
ing, the trajectory of ni(<) can become very complicated. 
While in this report we outline the general formalism, 
a detailed numerical study of the magnetization dynam- 
ics will be carried out elsewhere. In the remainder of 
this paper we discuss the critical current bias at which a 
collinear equilibrium configuration becomes unstable. 

Near a collinear configuration, Eq. |(SJ) simplifies to 
jcxch _ IPu™P/2. Let m (circularly) precesses around M 
with the FMR frequency w: m x dm/dt = wm x M x m. 
The total (projected) spin current in the Gilbert form 
then reads: 

m X I, X m = [%Ti/(87r) -|- lf''y{2ij)] m x ^ . (12) 
An instability is reached when the effective Gilbert damp- 



ing coefficient becomes negative. The critical bias is thus 
given by 

Is!!" - [.9^V(47r) + 2aoMsSd/{hj)] Tiw . (13) 

Neglecting the first term in the brackets of the above 
expression, we obtain result analogous to Slonczewski's^ 
while neglecting the second term, we get a condition sim- 
ilar to Berger's spin- wave emission criterionnSi The spin- 
pumping contribution (first term) is comparable with the 
intrinsic damping (second term) for films with thickness 
d of several nanometersj^iifi with the former dominating 
for very thin films. 

In summary, we have developed a general theoretical 
framework for the low-temperature magnetization dy- 
namics in small spin valves, unifying and extending pi- 
oneering work by SlonczewskiS and Berger^ The non- 
equilibrium spin torque induced by the bias current and 
the enhanced Gilbert constant due to the spin pump- 
ing must be treated on equal footing. When the mem- 
ory magnetic element is sufficiently thin (d < 10 nm), 
the nontrivial dependence of the damping on both the 
static and dynamic configurations of the system can im- 
portantly modify the magnetization dynamics. We de- 
rived the dependence of the Gilbert damping of on 
the precession-cone angle, which can also be measured by 
the FMRii^ Micromagnetic simulation codes should take 
these effects into account as the device and magnetic bit 
dimensions decrease down to the nanometer scale. 
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